NPS  ARCHIVE 
1969 
BOND,  W. 


INVESTIGATION  OF  DISTORTION  OF  DIVERS' 

SPEECH  USING  POWER  SPECTRAL  ESTIMATES 

BASED  ON  THE  FAST  FOURIER  TRANSFORM 


by 


Wi  1 1 iam  Howard  Bond 


DUDLEY  KNOX  LIBRARY 
MONTEREY,  CA  93943-5101 


United    States 
Naval   Postgraduate   School 


THESIS 


INVESTIGATION 

OF  DISTORTION  OF  DIVERS' 

SPEECH 

USING 

30WER  SPECTRAL  ESTIMATES 

BASED 

ON  THE 

FAST  FOURIER  TRANSFORM 

" 

by 

Will 

iam  Howard  Bond 
and 

James  Michael  Myatt 

June  1969 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  Is  unlimited. 


¥133498 


. 


LIBRARY 

KAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,.  CALIF.,  9394Q 


Thesis  by:   William  H.  Bond  and  James  M.  Myatt  entitled  "Investigation 
of  Distortion  of  Divers'  Speech  Using  Power  Spectral  Estimates  Based  on 
the  Fast  Fourier  Transform."  June  1969. 

ERRATA 

Page  Line  Change  To 

58  8  Fig.  23  a-c  Fig.  23 

61  6  Fig.  20  (a-c)  Fig.  22  a-c 

61     2nd  paragraph,  insert  as  first  line  --  "The  transient  spectrum 
with  the  oral  cavity  is  shown  in  Fig.  24." 

61     3rd  paragraph,  insert  as  first  line  --  "Fig.  25  shows  the  transient 
spectrum  at  80  ft.  for  the  first  2048  points." 


Investigation  of  Distortion  of  Divers' 

Speech  Using  Power  Spectral  Estimates 

Based  on  the  Fast  Fourier  Transform 


by 


William  Howard^Bond 

Major,  United  States  Marine  Corns 

B.S.,  University  of  Texas,  1957 

and 

James  Michael  Myatt 

Captain,  United  States  Marine  Corps 

B.S.,  Sam  Houston  State  College,  1963 


Submitted  in  partial  fulfillment  of  the 
reguirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
June  1969 


ABSTRACT 

The  nroblem  of  distortion  in  underwater  communications  peculiar 
to  free  divers  and  techniques  for  analysis  of  sneech  wave  forms  are 
discussed.  The  Fast  Fourier  Transform  alqorithm,  selected  to  analyze 
shifts  in  formant  frequencies  due  to  restricted  oral  cavities,  hinh 
ambient  pressures,  and  forced  speech  is  discussed.  The  Fast  Fourier 
Transform  is  used  to  analyze  a  vowel  sound  and  show  that  the  expected 
shifts  do  occur.  Recommendations  are  made  for  extendi nq  the  techniques 
to  all  non-noise  like  sounds  and  breathing  mixtures  other  than  compressed 
air. 
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I.   INTRODUCTION 

Man's  future  needs  dictate  the  development  of  the  vast  areas  of 
the  continental  shelf  for  food,  raw  materials,  and  Hvina  space.  To 
explore,  map,  and  colonize  this  unknown  a  myriad  of  new  equipment  must 
be  developed  for  explorers  of  this  "last  frontier."  Today  man  has  the 
capability  to  communicate  reliably  over  the  vast  reaches  of  soace;  but 
free  divers  (Divers  with  no  attachment  to  a  habitat  or  surface  station  - 
untethered)  operating  in  close  proximity  cannot  communicate  reliably  to 
each  other.  One  UDT  officer  summed  up  the  state  of  the  art  by  sayina, 
"reliable  communications  between  free  divers  consists  of  tugs  on  buddy 
lines."  Unsatisfactory  undersea  communications  result  from  poorly 
designed  oral  cavities,  high  ambient  pressures,  poor  electronics,  and 
exotic  gas  breathing  mixtures.  Current  equipment  has  been  further 
complicated  by  designers  unfamiliar  with  the  environment  and  hazards 
associated  with  the  free  diving  world;  few  designers  have  taken  into 
consideration  the  user's  emotional  response  to  the  apparatus  in  the 
underwater  environment. 

In  1959-60,  the  problem  of  inadequate  communications  between  free 
divers  was  identified  in  a  specific  operational  requirement  of  Navy  UDT 
and  Marine  Corps  Force  Reconnaissance  Divers.  Since  that  time,  the 
U.  S.  Navy  and  U.  S.  Marine  Corps  have  been  moving  on  separate  paths  to 
develop  a  reliable  communicator  for  the  free  diver.  Presently  available 
equipment  is  far  from  an  optimum  design  for  many  reasons  —  poor  overall 
intelligibility  being  of  prime  importance.  In  short,  the  state  of  the 
art  is  unacceptable.  This  unacceptable  condition,  personal  interest, 
and  the  requirement  of  the  Armed  Forces  for  more  sophisticated  devices 
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qenerated  the  oriqinal  tonic  of  this  research  -  an  investiaatlon  into 
sneech  intelligibility  of  free  diver  underwater  communicators.  After 
visitinn  with  nersons  enaaned  in  Navy  sponsored  work  in  this  realm,  it 
became  apparent  that  more  was  Involved  than  just  better  electronics. 
Inherent  in  improvina  underwater  communications  devices  1s  the  need 
for  better  understand!' nq  of  what  happens  to  speech  in  the  undersea 
environment.  This  fact  made  it  necessary  to  pursue  a  workinq  knowledqe 
of  the  disciplines  of  linquistlcs  and  speech  processinq.  Repeated 
conferences  with  the  Communications  Sciences  Laboratory,  University  of 
Florida  stimulated  the  curiosity  and  study  of  these  disciplines. 

Past  speech  processing  has  taken  the  form  of  qualitative  analysis 
usinq  analoq  equipment,  such  as  spectoqraphs,  and  intelligibility 
studies  usinq  human  listeners;  however,  few  quantitative  studies  have 
been  undertaken  to  isolate  the  major  distortion  problems.  Many  compar- 
ative studies  of  hard  line  and  free  diver  communication  systems  have 
been  and  are  beinq  carried  out,  primarily  usinq  human  speakers  and 
listeners  to  evaluate  system  intellioibility  usinq  controlled  word 
lists  [1].  The  spectoqraph  also  gives  a  measure  of  quantitative  anal- 
ysis to  highly  trained  speech  analysts. 

The  use  of  frequency  analysis  (a  common  tool  of  the  electrical 
enqineer)  for  simple  waveform  analysis  is  a  proven  technique.  In  the 
past,  this  technique,  transforming  the  siqnal  under  consideration  to 
the  frequency  domain  usinq  the  Fourier  transform  required  too  areat  an 
expenditure  of  time  and  human  effort  for  speech  waveforms.  The  advent 
of  diqital  computers  with  vast  memory  capabilities,  now  makes  this 
technique  a  feasible  method  for  speech  waveform  analysis.  Transfor- 
mation of  the  speech  waveforms  is  accomplished  usinq  the  discrete 
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Fourier  transform.  This  transformation  1s  reliable  if  the  siqnal  is 
band  limited  and  occurs  for  a  finite  period  of  time.  The  Fast  Fourier 
Transform  (FFT)  alqorithm  provides  a  rapid  method  for  calculatinq  the 
discrete  Fourier  transform  of  such  slqnals  [2].  This  alqorithm  has 
many  applications  to  dlqltal  filterlnq  and  control linq  real  time  proc- 
esses in  the  frequency  domain.  The  value  is  apparent  when  conslderinq 
siqnals  that  require  8,000  samples  to  provide  complete  coveraqe  of 
frequencies  of  interest.  It  provides  a  real  time  savlnq  of  300  to  1 
over  previous  alqorithms;  analysis  that  previously  took  hours  now  takes 
seconds. 

This  1s  the  approach  taken  1n  this  paper  -  1nvest1qation  of  dis- 
tortion in  small  underwater  communicators  usinq  the  FFT  to  process 
diver's  speech.  The  final  value  can  only  be  Imaqined;  the  analysis  can 
be  applied  to  isolate  distortion  in  helium  speech,  build  better  helium 
speech  unscramblers  and  help  to  prevent  loss  of  vital  undersea  explorers, 
such  as  Berry  Cannon. 

A.  FREE  DIVER  COMMUNICATION 

Whenever  man  qoes  underwater,  he  needs  communication.  Tasks  under- 
water are  made  easier  and  safer  1f  communications  exist  between  divers 
and/or  surface  stations.  The  military  diver  uses  communications  to 
carry  out  complicated  tasks,  to  penetrate  enemy  held  areas,  and  to  effect 
rendezvous  with  underwater  vehicles.  Reliable  underwater  communications 
do  not  exist  for  all  cases  of  interest  to  the  U.  S.  Navy  and  U.  S.  Marine 
Corps. 

Communications  between  free  divers  is  a  relatively  new  need;  but  as 
man  qoes  deeper  for  lonqer  periods  of  time  the  requirement  will  increase 
in  importance.  Sonar  has  been  used  for  communications  between 
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submersibles  and/or  surface  vessels  for  many  years;  small  sonar  sets 
are  now  performinq  the  same  service  between  divers,  submersibles  and/or 
surface  stations.  Many  problems  exist  in  the  diver's  case  that  are  not 
of  consideration  in  submarines.  The  more  important  include  the  addition 
of  restricted  oral  cavities  to  the  vocal  tract,  increased  ambient  pres- 
sures, increased  noise  levels,  and  the  uses  of  exotic  breathinq  mixtures 
1 .  Current  Equipment 

Initially,  efforts  to  provide  the  diver  with  communication 
involved  amDlifvinq  the  voice  and  projectinq  it  into  the  underwater 
environment.  Water  noises  near  transmission  frequencies,  low  power,  and 
exclusion  of  microphones  and  earphones  (for  low  cost)  contributed  to 
Door  intelligibility  and  short  ranqe  makinq  these  systems  impractical 
for  most  military  uses  [3],  Electromaqnetic  and  electric  field  trans- 
mission have  also  been  tried  with  limited  success  at  short  ranqes. 

Sonar  has  provided  the  most  successful  method  of  communication 
to  date.  Using  crystals  possessing  electrostrictive  properties,  a  sound 
pressure  or  acoustic  wave  is  generated  in  to  the  water  and  received  by  a 
similar  transducer.  Sets  of  this  design  operate  around  a  carrier  fre- 
quency of  8KHz-42KHz  (the  lower  carrier  producing  greater  ranqe  of 
transmission)  and  utilize  AM-SSB/SC  transmission  to  lower  noise  recep- 
tion and  increase  power.  The  Hydro-Products  811  and  Bendix  PflC-2  (see 
Figs.  1  and  2)  underwater  communicators  are  industry's  latest  efforts  in 
the  field.  Characteristics  of  these  sets  are  shown  in  Table  I. 

Chanqes  occur  in  these  sets  when  placed  on  the  diver.  The  omni- 
directional transducer  becomes  somewhat  directional  due  to  acoustical 
dampeninq  by  the  lungs,  tanks,  and  the  neoprene  wet  suit  worn  in  cold 
water.  Two  divers  operating  in  close  proximity  can  not  communicate  when 
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Hydro  Products  811 
Underwater  Communicator 

FIGURE  1 


II 


Bendix  PQC-2 
Underwater  Communicator 

FIGURE  2 


15 


TABLE  I 


SET 


HP  811 


POC-2 


Carrier  Frequency 

8.0875  KHz  * 
AM-SSB/SC 

8.0875  KHz  * 
AM-SSB/SC 

Electric  Power 
Into  Water 

.5  watt     * 

.6  watt    * 

Ranqe 

Nominal  800  meter  * 

500  meter  minimu 

Modes 

Voice 

Voice,  Tone,  and 
Interruoted  Tone 

Transducer 

Electrostrictive 
(barium  titamate) 

Electrostrictive 
(lead  zirconate) 

Pattern 

Omnidirectional 

Omnidirectional 

Microphone 

Oral  Cavity 
Microphone 

Bone  conduction 
(two) 

Receiver 

Bone  Conduction 

Bone  Conduction 

(one) 


(two) 


Intel  1 iqibil ity 


No  information 
(under  test  at 
U.  of  Fla.) 


No  information 
(under  test  at 
U.  of  Fla.) 


Bandwidth 


300-3000  Hz  * 
(2700  Hz) 


300-3000  Hz  * 
(2700  Hz) 


*  Manufacturer's  Data 
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the  lunq  cavity  is  1n  the  line  of  transducers.  This  is  a  serious 
defect  when  one  user  is  unaware  of  the  relative  location  of  the  second 
user. 

The  communicator  bandwidth  operates  on  the  telephone  principle  - 
300-3000  Hertz  bandwidth  will  provide  hiah  Intelligibility.  Studies 
indicate  for  shallow  depths  that  this  bandwidth  will  provide  intelli- 
gibility of  95+%  [1];  but  this  1s  not  necessarily  true  for  all  conditions 
encountered  by  the  diver. 

Evidence  points  to  bone  conduction  as  the  primary  source  of 
hearing  underwater  [4],  All  currently  available  communicators  utilize 
bone  conduction  reception;  one  of  the  tested  communicators  used  bone 
conduction  transmission  in  place  of  the  normal  microphone  in  the  oral 
cavity.  Some  tests  Indicate  a  definite  loss  1n  intelligibility  when 
using  bone  conduction  transmission;  but  no  quantitative  results  support 
these  conclusions.  The  use  of  a  single  hose  requlator  with  bone  con- 
duction transmission  caused  the  signal  to  be  obscured  by  bubble  noise 
passing  the  conductors. 

These  problems  are  not  unique  to  the  tested  communicators  but 
a  far  more  important  set  of  problems  are  those  common  to  all  communi- 
cation systems  for  free  divers. 
2.  Problem  Areas 

The  most  complex  problem  is  that  of  speech  distortion  caused  by 
the  addition  of  a  small  cavity  over  the  mouth.  To  date  little  basic 
research  has  been  attempted  in  an  effort  to  design  an  optimal  oral 
cavity.  The  design  must  balance  two  conflicting  variables  -  small 
volume  to  prevent  low  resonances  near  formants  or  articulated  sound  and 
buildup  of  carbon  dioxide,  and  large  volume  to  prevent  increased 
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pressure  from  affecting  the  speaker  in  articulation.  The  Bioenqionics ' 
Nautilus  oral  cavity  (see  Fiqure  3)  was  used  to  aid  speaker  articula- 
tion in  this  research.  This  oral  cavity  has  been  used  successfully  to 
increase  intelligibility  of  systems  in  other  tests  and  nrovides  an 
excellent  trade  off  in  cavity  desiqn  and  diver  safety  [5]. 

Increased  operating  depth  results  in  increases  of  ambient  pres- 
sures .  This  chanqes  the  characteristics  of  the  resonant  cavities  and 
the  density  of  the  breathing  qas.  When  a  diver  is  ODeratinq  at  pres- 
sures of  four  atmospheres  (100  feet)  or  qreater,  the  voice  attains  a 
nasal  quality  and  causes  considerable  loss  of  intelligibility  [5]. 

When  exotic  breathinq  qases  are  substituted  for  the  primarily 
nitroqen-oxyqen  qas,  the  voice  attains  a  quality  best  described  as  the 
"Donald  Duck"  effect.  The  use  of  exotic  gases  causes  formant  shifts 
due  to  the  increased  sound  velocities.  Unfortunately  the  shifts  are 
non-linear;  this  makes  correction  a  complex  Droblem.  The  amount  of 
distortion  due  to  these  effects  is  under  study  by  conventional  tech- 
niques by  the  Communications  Sciences  Laboratory,  University  of  Florida 
[1],  and  other  aqencies  [6,7,8]. 

Physioloqical  problems  are  of  interest  to  ensure  that  human 
enqineerinq  of  communicators  are  acceptable  to  the  diver  when  ODeratinq 
in  an  alien  environment  and  commensurate  with  desired  safety.  The 
diver  feels  an  almost  personal  relationship  with  his  equipment.  The 
communicator  must  require  little  or  no  chanqe  in  equipment  for  the  diver 
to  accept  it  unhesitatingly.  A  prime  example  is  the  addition  of  the 
oral  cavity,  a  much  needed  item  for  proper  articulation.  The  original 
desiqn  of  the  Bendix  POC-2  excluded  use  of  a  cavity  for  this  reason. 
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Desinn  of  workable  underwater  communicators  is  feasible  ;  but  a 
qreat  deal  of  research  is  necessary  before  a  sinqle  system  will  over- 
come the  vast  number  of  problems.  The  most  critical  item  1s  development 
of  a  method  to  solve  the  deep  depth  and  exotic  qas  speech  distortion  and 
then  microminaturlzation  of  the  "unscrambler"  for  Inclusion  in  the  free 
diver's  equipment. 

B.   SPEECH  PROCESSING 

To  the  engineer  engaqed  in  speech  research,  "speech  processinq" 
implies  synthesis  and  analysis,  automatic  recognition  and  speech  com- 
pression. The  linguist  will  speak  of  the  acoustical  model  in  speech 
production  or  perception  and/or  various  articulation  and  intelligibil- 
ity testing.  The  language  of  speech  research  is  foreign  to  the  beqinner; 
therefore,  Appendix  A  Includes  a  glossary  of  speech  and  linguistic  terms. 

As  previously  mentioned,  the  requirement  exists  for  a  technigue 
that  will  enable  the  speech  analyst  to  accurately  determine  the  effects 
of  an  alien  environment,  such  as  that  encountered  by  the  diver,  on  the 
human  voice. 

Some  of  the  technigues  used  for  speech  analysis  are  summarized  with 
a  brief  description  of  the  new  "wrinkle"  in  analysis  following. 

1 .  Analog  Techniques 

a.  Recognition  Analysis 

The  vowel  sounds  have  been  the  subject  of  intense  research. 
The  most  common  result  has  been  the  classification  of  the  vowels  by 
their  formants,  primarily  the  relation  of  the  formants  to  each  other, 
and  has  been  labeled  "Visible  Speech"  [9].  The  device  most  commonly 
used  is  the  "acoustic  spectogram"  where  a  visual  representation  of  the 
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spoken  messaqe  or  phoneme  is  displayed  in  freauency  versus  time. 
The  intensity,  a  third  parameter,  1s  displayed  by  the  darkness  of  the 
plot.  Thus,  the  transient  spectrum  can  be  viewed.  A  spectooram  of  the 
word  "at"  1s  shown  in  Fiq.  4.  The  chief  disadvantaae  of  the  spectoaram 
is  the  lack  of  resolution. 

b.  Power  Spectrum  Measurements 

The  analoq  spectrum  analyzer  is  a  convenient  tool  to  view 
the  spectrum  of  a  signal  instantaneously.  It  works  on  the  principle 
of  a  constant  width  bandpass  filter  with  a  variable  center  freauency; 
however,  hiqh  precision  measurements  are  not  possible.  Another  tech- 
nique that  borders  on  beinq  classed  as  diqital  is  the  use  of  a  bank  of 
narrow  bandpass  filters  with  increasinq  bandwidth  as  their  center  fre- 
quency increases.  Until  recently,  size  and  cost  have  limited  their 
use,  with  resolution  still  the  determininq  factor. 

c.  Waveform  Analysis 

The  fact  that  vowel  sounds  can  be  characterized  by  their 
waveform  asymmetry  has  been  noted.  Pairs  of  phonemes  can  be  compared 
and  separated  on  the  basis  of  the  magnitude  and  polarity  of  the  output 
of  a  demodulator  [10]. 

d.  Articulation  and  Intelligibility  Tests 

To  evaluate  communication  systems,  the  articulation  tests 
were  widely  used  during  World  War  II.  The  results  of  the  tests  aided 
hardware  design.  Intelligibility  testing  takes  the  form  of  readinq 
from  phonetically  balanced  word  lists  and  assigninq  a  percentaqe  score 
for  the  correct  responses  of  human  listeners.  The  consistency  of  such 
tests  demonstrates  their  reliability  and  makes  them  the  most  commonly 
used;  however,  the  search  for  optimum  lists  of  sounds,  words  and  phrases 
continues  [5]. 

21 


-a 

s- 
o 


+-> 

4- 
O 


s- 
cr. 
o 
+J 
o 
a> 

Q- 


o 

o 

o 

O 

O 

CD 

O 

o 

o 

o 

l£> 

CM 

CO 

J" 

O 

CD 

CM 

CO 

J- 

i—* 

en 

l£> 

J" 

CM 

en 

r~- 

J- 

CM 

CM 

i-H 

■H 

■— > 

»— < 

(zh)  ADNanOa^d 


22 


2.  Digital  Techniques 

The  analysis  of  speech  by  digital  techniques  owes  a  great 
deal  of  its  popularity  to  the  interest  in  the  development  of  the 
"vocoder",  or  voice  coder.  Linguists  and  psycholinguists  have  been 
a  trifle  reluctant  to  master  the  techniques  of  the  digital  world  [11]. 

Digital  Spectrum  analyzers  prior  to  1965  consisted  of  narrow 
band  digital  filters.  The  number  of  digital  filters  necessary  for 
sufficient  bandwidth  and  resolution  required  extensive  memory  and 
computing  time;  an  excessive  amount  for  speech  waveforms. 

a.  The  Fast  Fourier  Transform 

The  caDability  for  high  resolution  spectral  analysis  was 
not  available  until  recently  discovered  algorithms  were  imnlemented 
[12].  The  algorithm  of  interest  in  this  paper  is  the  Fast  Fourier 
Transform  (FFT)  [2].   It  is  an  efficient  algorithm  for  computing 

the  "discrete  Fourier  transform".  The  number  of  operations  is  reduced 

2 
from  N  to  2Nlog2N,  where  N  is  the  number  of  samples.  Thus,  it  can 

be  used  for  offline  spectral  analysis  of  processes.  The  advantage 

of  the  FFT  is  illustrated  by  comparing  the  number  of  operations  re- 

2 
quired  for  N  =  8192  samples.  N  is  in  excess  of  67  million  compared 

to  213  thousand  for  the  FFT.  Such  a  saving  in  ooerations  vastly 

reduces  the  computing  time. 

b.  Spectral  Measurements 

With  the  FFT  as  a  tool  available  in  generating  the  Fourier 
coefficients,  several  avenues  of  analysis  are  plausible.  The  avenue 
chosen  is  the  power  spectrum.  Other  spectra,  such  as  the  cesntrum, 
are  available  but  were  not  considered.  The  FFT  lends  itself  to  such 
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uses  and  could  be  apnlied  in  future  research  in  this  area.  Bv  claiminn 
that  the  power  spectrum  is  available,  it  must  be  clear  that  the  spec- 
trum is  only  an  estimate.  How  good  the  estimate  is  depends  upon  the 
methods  of  implementing  the  alqorithm.  First,  the  siqnal  is  of  finite 
lenqth,  hence  the  number  of  samnles  is  finite.  Second,  the  siqnal 
is  bandlimited.  The  assumption  is  made  that  the  siqnal  is  a  member 
of  a  stationary  process. 
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II.   DATA  COLLECTION 

After  consultation  with  authorities  in  sneech  nrocessinq  and 
underwater  communications,  and  after  considering  time,  funds,  and 
human  effort  available,  the  decision  was  made  to  proceed  to  collect 
data  to  suDport  two  independent  projects.  This  decision  was  influ- 
enced by  the  loan,  free  of  cost,  of  the  two  underwater  communicators 
used  in  the  data  collection.  Certain  control  sounds  and  words,  shown 
in  Table  II, were  recorded  to  support  the  digital  speech  processing. 
Two  vowel  and  four  consonants  were  arbitrarily  chosen.  One  sound 
of  each  main  subgroun  was  included  [9].  This  data  was  orininallv 
intended  to  assist  in  isolating  distortion  in  the  two  communicators. 
The  second  project  was  designed  to  assist  in  judgina  the  relative 
intelligibility  of  the  two  systems  under  various  conditions.  Two 
Campbell  word  lists  were  recorded  by  each  diver.  The  Camobell  word 
lists  are  monosyballic,  Dhonetically  balanced  lists  of  25  words  each, 
selected  randomly  from  a  master  list  of  500  words.  Table  III  is  an 
example  of  such  a  list. 

Oriqinal  nlans  called  for  data  to  be  taken  in  two  locations,  an 
anechoic  pool  free  of  outside  interference  and  the  ocean  to  nrovide 
a  noisy  environment  with  varying  conditions  of  range  and  depth.  Sea 
conditions  and  boat  availability  in  January  and  February,  1969 
allowed  only  eight  (8)  days  to  be  used  for  takinci  data;  althouqh  the 
effort  was  made  daily.  Finally  after  30%  of  the  desired  data  had  been 
collected,  theft  of  equipment,  lack  of  a  stable  platform  in  heavy  seas, 
poor  ocean  visibility,  and  shortage  of  personnel  dictated  a  change  to 
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Phoneme 


TABLE  II 
Control  Sounds  and  Words 

Word 


Class  of  Phoneme 


/AE/ 

AT 

Vowel 

Itl 

MET 

Vowel 

/?/ 

PAY 

Voiceless  Ston 

/G/ 

GO 

Voiced  Ston 

/// 

SHOULD 

Voiceless  Fricative 

/F/ 

FIRST 

Voiceless  Fricative 

/TH/ 

THIN 

Voiceless  Fricative 

/v/ 

VOTE 

Voiced  Fricative 

TABLE  III 
Sample  Campbell  List  (P-7) 


DOLLS 

NUTS 

BOOK 

JUMP 

ART 

CARVE 

OR 

ARE 


AIM 

COOK 

CHAIR 

AID 

SMOOTH 

AM 

ARM 

ACE 

POOR 

UP 

HIGH 

YES 

SKIN 

THREE 

LIE 

HURT 

GIVE 
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a  quiet,  fresh  water  environment  havinq  a  stable  platform.  Two  under- 
water communicators,  the  Hydro-Products  811  and  the  Bendix  POC-2  were 
used  to  transmit  the  data  to  an  AN/WfiC-lA  communication  set.  Charact- 
eristics of  this  set  appear  in  Table  IV.  This  receiver  fed  a  Precision 
Instruments  DA  6200  tape  recorder.  The  data  was  recorded  at  37.5  ios 
and  stored  on  maqnetic  tape.  All  recordinq  used  frequency  modulation 
to  allow  for  subsequent  linear  frequency  translation  by  chanqinq  tane 
recorder  sneed. 

A.  ANECHOIC  POOL 

The  anechoic  pool  located  in  Room  025  of  Snanaqel  Hall,  Naval 
Postqraduate  School,  was  utilized  to  provide  the  noise  free  environ- 
ment. The  diver  was  dressed  in  normal  divinq  equipment  as  shown  in 
Fig.  5  and  6.  Initially  the  diver  sooke  into  free  soace  as  the 
control  sounds,  control  words,  and  two  Campbell  word  lists  were 
recorded.  Then  an  oral  cavity  was  placed  in  position  and  the  proce- 
dure repeated.  A  block  diaqram  of  the  test  setup  is  shown  in  Fiq.  7. 

B.  OCEAN  RANGE 

The  ocean  ranqe  is  located  in  Monterey  Bay  1500  meters  off  Del  Monte 
Beach.  The  location  is  shown  in  Fiq.  8.  The  underwater  portion  of 

this  ranqe  is  still  located  as  indicated.  The  ranqe  is  located  in  100 

feet  of  water  over  a  flat  sandy  bottom  to  minimize  bottom  reflection 

interference.  Fiq.  9  and  10  show  the  ranqe  layout.  A  diver  platform 

is  located  at  the  main  buoy  to  allow  orientation  of  the  diver  toward 

the  receivinq  transducer.  This  orientinq  device  orevents  siqnal  loss 

due  to  attenuation  by  the  lunq  cavities,  wet  suits,  and  other  items  of 

diver  dress.  Fiqure  11  shows  a  diver  in  the  platform  ready  to 
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TABLE  IV 
Characteristics  of  AN/WQC-1 A  Communication  Set 


OPERATION 
FREQUENCY 
MODES 
TRANSDUCER 

BANDWIDTH 


AM/SSB-SC 

8.3-11   KHz 

Receive,  Transmit,  Tone 

Omnidirectional  w/o  Baffle 

15  db  down  60°-330°  w/  Baffle 

300-3000  Hz 
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Diver  in  Anechoic  Pool 
Without  Oral  Cavity 

FIGURE  5 


Divsr  in  Anechoic  Pool 
With  Oral  Cavity 

FIGURE  6 


on 


O         uu 
O         Q 

cvj       a: 
vo  uj  o 

Q.  h-  Q£ 


I 


I        UJ 

Cr    uj 

3     Ct 


I 


Q£ 

UJ 

C_) 

ID 

<c 

Q 

t— 

OO 

1 

Z 

(_> 

<x. 

cr 

or 

3 

t— 

Q. 
■»-> 

W>    O 

O)   Q. 


o  o 

E    O 

<T5    Ol 

s-  c 
cn«=c 


OC 
C3 


o 

CO 


30 


o 

•r— 
-t-> 
03 
O 

o 


a» 
en 

c 

03 
03 

at 
o 
o 


CL 


CO 


CD 


31 


o- 


* 


u 


o 

o 


CO 


_  J 


Q> 


O 


R 


1= 


J 


O    li 


o 


l\ 


o 

LO 

s   i 


1   \ 


2 
cu 


Q. 
O 


e 


0) 

o 

o 


C7» 


o 


o 

o 
o 


co 


CO 


o 


32 


<r 


■ 


<■ 


CC 

IS) 

0 

\- 

ZT 

< 

0 

o 

C_) 

_J 

^ < 

u_ 

1— 

1- 

oo 

z: 

<c 

•a. 

_J 

>- 

-J 

0 

<r 

ID 

CO 

CO 

=*fc 

=*= 

C\J 

0 

VO 

IT) 

^— 

<  o 


CO 


W) 


4-> 


<u 

> 


A 


E 
O 


L 


OJ 


rr~j 


<■ 


r 


1 


<]   ™     1 


<J- 


-5'      ^     [__1_J 

^20'  >i    30'   -4-<30' 
K  100' 


A 

O 

/  • 

O 

in 

O) 

> 

O 

E 

0) 

t— 

O 

■a 

O 

•r— 

LU 

C\J 

U~> 

cc 

Z3 

fl 

O 

O) 

1— 1 

CD 

U_ 

c 

10 

Q£ 

C 

lO 

(1) 

O 

O 

\ 


o 
o 


5V  \ 


33 


Diver  in  Orientation  Platform 
FIGURE  11 
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communicate.  Finure  12  is  a  block  diaqram  of  the  test  set-up.  Note 
that  the  tane  recorder  is  mounted  on  an  unstable  nlatform  (the  boat) 
and  must  remain  within  short  ranqe  (about  15  meters)  of  the  desired 
buov  ranqe  marker.  Table  V  indicates  the  desired  conditions  for 
recordinq  data. 

C.   FRESH  WATER  RANGE 

The  fresh  water  ranqe  was  set  up  in  San  Vincente  Reservoir,  San 
Dieqo  County,  California.  Fiqure  13  shows  the  site.   It  differs  from 
the  previously  described  ocean  ranqe  in  four  ways. 

1.  The  granite  walls  and  concrete  dam  face  made  receivinq 
transducer  orientation  critical  to  prevent  side  reflections  and 
minimize  reverberation.  This  was  accomplished  by  susnendina  the 
transducer  on  one  line  and  usinq  two  quidelines  to  insure  exact 
posi tioninq. 

2.  The  tape  recorder  was  olaced  on  a  stable  barqe  and  the 
diver  moved  to  different  locations,  thus  reversinq  the  roles  in  the 
ocean. 

3.  Ouiet,  fresh  water  replaced  the  noisy  ocean  environment. 

4.  Diver  orientation  was  accomnlished  usinn  a  comoass. 
Table  V  shows  the  desired  conditions  for  data  collection.  Eiqht 
days  were  required  for  movement  to  the  reservoir,  set  un  on  the 
barqe,  collection  of  data  and  return  to  Monterey,  California. 
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San  Vicente  Reservoir  Test  Site 


FIGURE  13 
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TABLE  V 


Recordinq  Conditions  for  Ocean  and  Fresh  Water  Ranges 


RANGE 

10 
10 
10 

100 
100 
100 

300 
500 


DIVER 

BOND 

MYATT 

SET 

P0C-2 

HP  811 

PQC-2 

HP 

DEPTH 

20 

X  Y 

Y 

X  Y 

Y 

50 

X  Y 

Y 

X  Y 

Y 

80 

X  Y 

Y 

X  Y 

Y 

20 

X  Y 

Y 

X  Y 

Y 

50 

X  Y 

Y 

X  Y 

Y 

80 

X  Y 

Y 

X  Y 

Y 

50 

Y 

Y 

Y 

Y 

50 

Y 

Y 

Y 

Y 

X  -  DATA  RECORDED  IN  OCEAN  RANGE 

Y  -  DATA  RECORDED  IN  FRESH-WATER  RANGE 
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III.  DATA  PROCESSING 

The  instruments  available  at  the  Naval  Postgraduate  School  for 
the  orocessinq  of  the  data  include  the  IBM  360/67  comouter,  the  SDS 
9300  medium  sized  comDuter  with  a  hybrid  interface  to  a  COMCOR  5000 
analoq  comouter. 

A.  ANALOG  TO  DIGITAL  CONVERSION 

The  method  of  Drocessing  the  data  deoended  largelv  on  the  capa- 
bilities of  the  hybrid  comouter  setup  of  the  COMCOR  5000  and  the  SDS 
9300.   In  order  to  digitize  the  data  from  the  tapes,  the  sanmlinn  rate 
of  the  analoq  to  digital  (A-D)  converter  had  to  be  at  least  twice  as 
high  as  the  largest  frequency  of  interest  (the  Nyquist  frequency). 
Unfortunately,  the  present  write  scheme  of  the  SDS  9300  limits  the 
samnle  rate  to  2500  Hertz.  This  virtually  eliminates  the  use  of  the 
A-D  converter  for  the  entire  audioband.   It  was  determined  that  the 
frequencies  of  interest  for  this  exoeriment  sliqhtly  exceed  the  nass- 
bands  of  the  two  communicators,  i.e.,  300-3000  hertz;  a  passband  of 
200-4000  hertz  was  used. 

The  samplinq  frequency  dilemna  was  circumvented  bv  utilizing  the 
features  of  the  FM  tape  recorder.  The  PI  6200  instrumentation  tape 
recorder  will  record  at  three  sneeds  -  37.5  ins,  3.75  ips,  and  .375 
i'ds.  The  center  frequency  of  the  carrier  is  qoverned  by  the  choice  of 
the  tane  sneed  -  50,000  hertz,  5,000  hertz  and  500  hertz,  respectively . 
Suitable  lownass  filters  are  switched  into  the  circuitry  for  these 
speeds.  At  37.5  ips  the  tape  recorder  input  yields  0  db  qain  from 
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dc  to  10,000  hertz,  with  a  signal  to  noise  ratio  of  42  db.  By  olavina 
back  at  3.75  i'ds,  the  discriminator  frequency  is  divided  by  a  factor 
ten,  as  are  the  siqnal  frequencies.  This  linear  frequencv  translation 
allows  a  samnle  rate  of  800  hertz  to  cover  the  translated  audioband. 

The  choice  of  parameters  for  sampling  the  data  was  hampered  by 
the  necessity  for  trade-off  in  frequency  and  time  resolution,  often 
referred  to  as  the  analogy  to  Heisenberqs  "uncertaintv  nrinciolc"  [13] 
Innorina  statistical  stability  and  those  terms  discussed  in  references 
on  nower  spectra,  the  frequency  resolution,  Af,  is  defined  as: 

Af  =  fs/N 
where  f   =  sampling  frequency 

N   =  number  of  samples 
The  dilemna  is  obvious.  The  larqer  the  number  of  samples,  the  smaller 
the  separation  between  spectral  lines,  hence,  qreater  resolution.  The 
larner  N  also  extends  the  time  duration  the  sinnal  must  be  samnled, 
hence  a  loss  of  time  resolution.  To  prevent  aliasinq,  the  samnlinq 
frequency  must  be  at  least  twice  the  highest  frequency  of  the  siqnal. 
The  larger  samolinq  frequency  reduces  the  frequencv  resolution. 

The  solution  was  to  investigate  the  data  for  both  time  and  fre- 
quency resolution.  The  time  resolution  was  increased  by  using  a  lower 
number  of  samples  for  calculating  an  "evolutionary  spectrum"  and  a 
larger  number  of  samples  for  a  spectrum  of  high  resolution. 

The  samnle  rate  was  selected  to  be  2043  hertz.  This  frequency  was 
selected  because  the  FFT  alqorithm  is  most  efficient  when  f  is  some 
number  that  is  a  power  of  two.   It  is  more  than  five  times  the  hiqhest 
frequencv  of  interest;  i.e.,  400  hertz  (4,000  hertz  translated  by  a 
factor  of  ten) . 
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The  number  of  samples  taken  must  necessarily  be  a  number  that  is 
a  power  of  two  for  the  FFT  algorithm.  Thus,  the  number  of  samnles  taken 
for  each  record  was  8192.  This  allowed  the  siqnal  to  be  sampled  for 
a  period  of  0.4  seconds.  The  maximum  resolution  possible  is  0.25 
hertz,  translated  by  the  factor  of  ten  to  2.5  hertz,  real  time. 

To  samnle  the  data,  the  tape  recorder  output  was  filtered  by  a 
Krohn-hite  adjustable  bandpass  filter,  Model  310  A.  The  character- 
istics of  this  filter  are  shown  in  Table  VI.  The  frequency  limits 
were  20  and  400  hertz.  The  siqnal  was  then  amplified  by  a  factor  of 
30,  and  diqitized  with  the  clock  rate  controlled  bv  the  Wavetek  116 
oscillator.  The  block  diagram  of  Fiq.  14  shows  the  samplinq  setun. 

The  siqnal  was  monitored  by  audio  means  and  visuallv  on  an  oscillo- 
scope to  determine  the  proper  time  to  commence  digitizing.  Herein 
lies  a  fundamental  weakness  in  this  experiment;  it  will  be  discussed 
more  fully  in  the  analysis  section.  Briefly,  no  control  was  used  to 
ensure  that  the  sampling  process  started  at  the  same  instant  the 
phoneme  commenced. 

A  vast  amount  of  data  was  recorded  for  analysis.  With  time 
becominq  a  factor,  it  was  necessary  to  limit  the  number  of  phonemes 
examined. 

The  phoneme  selected  was  |AE|,  the  vowel  sound  in  "at"  -  a  front 
vowel.  The  same  phoneme  was  sampled  for  the  followinq  conditions  - 
each  communicator  at  a  ranqe  of  ten  meters  and  depth  of  80  feet,  and 
each  communicator  in  the  anechoic  pool  with  and  without  oral  cavity. 

The  sampled  data  from  the  SDS  9300  was  placed  on  maqnetic  tape. 
The  output  was  seven  track  in  24  bit  octal  format.  The  samnle  routine 
was  an  in-house  routine  developed  by  former  NPGS  students  [14,15]. 
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TABLE  VI 
Characteristics  of  Model  310A  Krohn-hite  Filter 


1 

0 

-1 
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KROHN-HITE  MODEL  31 OC 
NORMALIZED  FILTER  RESPONSE 

Attenuation:  24  db/octave  below  -12db  point 

Input  Characteristics:  6  Megohms  in  parallel  with  50  ppf 

5  volts  rms  maximum 

Output  Characteristics:  200  ohms,  12  milliwatts 
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B.  INTERFACING 
The  output  of  the  SDS  9300  is  written  in  seven  track.  As  is  the 

case  with  all  comouters  when  tryinq  to  mate  with  the  IBM  360,  the 
interface  is  not  as  smoothly  refined  as  miqht  be  exnected.  The  seven 
track  tape  must  be  converted  to  nine  tracks  and  the  24  bit,  octal 
format  converted  to  32  bit,  hexadecimal  format.  The  oroqram  required 
is  1 isted  in  Appendix  B. 

The  tape  drives  of  the  IBM  360  required  perseverance  on  the  part 
of  the  onerators .  The  tane  from  the  SDS  9300  was  accented  rouqhly 
one  out  of  five  tries. 

No  attempt  was  made  to  store  this  data  on  a  direct  access  device. 
However, this  seems  desirable  to  conserve  computer  time. 

C.  DIGITAL  PROCESSING 
The  data  was  stored  on  maqnetic  tape  in  32  bit  words.  The  sub- 
routine Harm,  part  of  the  Standard  Scientific  Packaqe  for  the  IBM  360/67, 
was  used  to  calculate  the  estimates  for  the  Fourier  coefficients  of 

the  siqnal .  This  routine  is  based  on  the  Cooley-Tukey  alqorithm. 

As  previously  mentioned,  the  entire  8192  data  noints  could  be  used 
to  compute  a  soectrum  of  high  frequency  resolution  but  with  time  uncer- 
tainties. Alternately,  only  a  portion  of  the  data  points  could  be 
used  with  decreased  frequency  resolution  but  with  increased  certainty 
in  the  time  domain. 

The  techniques  complement  each  other.  The  technique  chosen  comnuted 
a)  a  power  spectrum  using  the  entire  8192  points,  and,  b)  to  assist  in  th 
analysis  of  this  spectrum, computed  a  transient  spectrum  for  the  first 
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2048  points.  The  two  spectra  are  titled  "spectrum"  and  "transient 
snectrum."  The  snectrum  has  a  resolution  of  2.5  hertz.  The  transient 
spectrum  has  a  frequency  resolution  of  ten  hertz. 

The  ontion  exists  whether  or  not  a  smoothinq  function  is  necessary. 
The  need  for  a  smoothinq  function  miqht  be  described  by  mental lv 
visual izinq  the  side  lobes  of  the  sample  function.  These  side  lobes 
cause  "leakaqe"  from  one  sample  to  another.  A  smoothinq  function 

such  as  the  now  familiar  "Hanninq  spectral  window"  reduces  the  side 

-1         -3 
lobe  leakaqe  from  (w-w1)   to  (w-w1)   [16,17].  Whether  or  not  the 

smoothinq  function  is  desirable  depends  on  the  tyDe  of  siqnal  .  The 
literature  says  that  noise-like  siqnals  with  relatively  flat  spectrum 
are  not  suitable  for  smoothinq  functions,  but  when  the  siqnal  is  non- 
noise-like  and  the  nosition  of  siqnal  enerqv  in  the  frequency  domain 
is  more  imDortant  than  the  amount  of  enerqy,  then  a  smoothinq  function 
is  anDroDriate.  The  "Hanninq  window"  was  used.  The  estimates  of  the 
Fourier  coefficients  were  oDerated  on  as  follows  (for  N  coefficients), 

A(k)  =  -3*a(k+l)  -'4  a(k-l)  +  \  a(k) 

B(k)  =  -h  b(k+l)  -\    b(k-l)  +  \  b(k) 

A(0)  =   %(a(0)  -  a(l)) 

B(0)  =0  k  =  1,2,3,4 

At  the  exnense  of  a  reduction  in  "statistical  stability",  this 
window  was  used;  its  time  domain  equivalent  beinq, 

X(j)  =  \  (l-cos(2  j/n))  x(j)    j  =  0,1,2 n-1 

A  second  smoothinq  oneration  was  performed  after  the  coefficients  of 
the  power  spectrum  were  computed. 

P(k)2  =  A(k)2  +  B(k)2     k  =  0,1,2 N/2 

P(k)  =   ;4(p(k+l)  +  n(k-l))  +  \n(k)      k  =  0,1,2 N/2 
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This  is  aqain  the  familiar  "Hanninq  window"  [17].  Since  the  main 
concern  is  that  the  estimates  are  anpronriate  to  the  nominal  fre- 
quencies, the  linearly  hanned  spectrum  suits  its  purnose. 

The  Dower  spectra  were  nlotted  on  a  loqarithmic  scale,  allowinq 
a  forced  scale  suitable  for  the  entire  dynamic  rancie  of  data  noints 
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IV.  ANALYSIS 

The  phoneme  selected  for  study  was  not  an  arbitrary  choice.  The 
noise-like  sounds  such  as  fricatives  or  stops  do  not  lend  themselves 
to  Fourier  analysis  [18].  The  vowel  sounds,  however,  have  a  steady 
state  quality,  with  slowlv  chanqinq  transient  SDectra.  The  front 
vowel,  /AE/  as  in  "at",  demonstrates  a  well  defined  snectrum  with 
three  important  forma nt  frequencies.  The  literature  shows  qeneral 
agreement  on  the  values  of  these  formants : 

F,  =  660  hertz 

F2  =  1720  hertz 

F3  =  2410  hertz  [19] 
Associated  with  these  formants  are  bandwidths  of  60,  100,  and  120 
hertz,  respectively.  One  reference  shows  the  formants  to  be  somewhat 
lower. 

F,  =  579  hertz 

F2  =  1691  hertz  [20] 

A.   FORCED  AND  NORMAL  SPEECH 

In  order  to  establish  some  sort  of  reference,  a  diver's  voice  was 
recorded  on  a  Kay  Missilyzer  soectoqranh.  Fiqure  15  shows  the  front 
vowel  /AE/  recorded  in  an  anechoic  chamber.  An  Altec  21BR150-2  micro- 
phone was  used  and  the  diver  spoke  in  a  normal  voice.  Fiqure  16  is 
the  same  sound  spoken  in  a  loud,  forceful  manner,  just  as  if  the  diver 
were  using  an  underwater  communicator.  The  requirement  for  the  forced 
sneech  to  increase  intelligibility  has  been  documented  [1]. 
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The  formant  frequencies  for  the  normal  sneech  are  shown  with 
arrows  in  Fiqure  15.  They  generally  aqree  with  the  values  cited 
previously.  Fiqure  16  shows  the  formants  greatly  increased.  The 
effort  to  increase  the  level  of  sound  has  caused  the  diver  to  lift  the 
pitch  of  his  voice  in  an  unnatural  manner.  The  soectrum  is  not  constant 
but  rises  raDi'dly  to  a  maximum  value  for  each  formant  and  then  slowlv 
decays.  This  chanqe  in  the  spectrum  is  of  major  importance.  Further 
verification  of  this  difference  in  the  forced  and  normal  speech  cower 
spectral  estimates  is  borne  out  by  Fiqures  17  (a-c)  and  18  (a-c). 
The  normal  speech  power  spectral  outout  closely  approximates  a  line 
spectra  where  shifts  can  be  easily  measured. 

B.  FREQUENCY  TRANSLATION 

To  check  the  frequency  translation  of  the  PI -6200  tane  recorder, 
a  384  hertz  tone  was  recorded  at  37.5  ins  and  samnled  at  the  rate 
of  2048  hertz.  The  power  spectrum  is  shown  in  Fiqure  19.  Fiqure  20 
shows  the  power  spectrum  when  the  tape  recorder  speed  is  slowed  to 
3.75  ips.  The  expected  chanqe  in  frequency  by  a  factor  of  ten  occurred. 

C.  POWER  SPECTRAL  ESTIMATES 
1 .  Without  Oral  Cavity 

A  sample  spectrum  for  the  recordinq  made  usinq  the  commun- 
icators without  oral  cavitv  is  shown  in  Fiqure  21 (a-c).  The  spectrum  is 
cluttered  with  power  at  frequencies  in  the  neighborhood  of  the  for- 
mants. This  clutter  is  sufficient  to  prevent  accurate  measurement  of 
the  formants;  but  should  be  anticipated  in  liqht  of  the  observations 
made  concerning  the  rapid  changes  in  the  soectrum  in  Figure  16.  This 
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FIGURE  17(a) 
Normal  Speech  Spectrum  of  /AE/,  0-1000  Hz. 
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FIGURE  18(a) 
Forced  Speech  Spectrum  of  /AE/,  0-1000  Hz. 
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FIGURE  17(b) 
Normal  Speech  Spectrum  of  /AE/,  1-2  KHz. 
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FIGURE  18(b) 
Forced  Speech  Spectrum  of  /AE/,  1-2  KHz. 
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FIGURE  17(c) 
Normal  Speech  Spectrum  of  /AE/,  2-3  KHz. 
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FIGURE  18(c) 
Forced  Speech  Spectrum  of  /AE/,  2-3  KHz. 
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FIGURE  19 
384  Hz.  Tone  at  37.5  ips 
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FIGURE  20 
384  Hz.  Tone  at  3.75  ips 
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FIGURE  21(a) 

Phoneme  /AE/,  Without  Oral  Cavity,  0-1  KHz. 
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FIGURE  22(a) 
Phoneme/AE/,  With  Oral  Cavity,  0-1  KHz. 
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FIG.  21(b) 
Phoneme  /AE/,  Without  Oral  Cavity,  1-2  KHz 
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FIG.  22(b) 
Phoneme  /AE/,  With  Oral  Cavity,  1-2  KHz. 
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FIG.  21(c) 
Phoneme  Without  Oral  Cavity,  2-3  KHz 
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FIG.  22(c) 
Phoneme  /AE/,  With  Oral  Cavity,  2-3  KHz. 
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snectrum  was  estimated  using  8192  samole  points.  Theoretically,  the 
frequency  resolution  is  2.5  hertz.  In  reality,  only  a  rough  estimate 
can  be  made  to  determine  the  formants  within  50  hertz.  Blaminq  the 
poor  results  on  the  manner  in  which  the  diver  forced  the  sound,  the 
next  step  was  to  compute  an  'evolutionary  spectrum',  (transient 
spectrum).  Using  the  first  2048  samnle  points,  the  theoretical  resolu- 
tion is  ten  hertz.  The  transient  spectrum  for  the  first  2048  points 
is  shown  in  Figure  23a-c.  This  covers  a  signal  for  0.1  second.  The 
obvious  blunder  is  that  no  technique  was  used  to  ensure  that  the 
sampling  process  beqan  precisely  at  the  beginning  of  the  sound.  This 
is  critical  in  the  transient  spectrum,  in  that  an  additional  uncer- 
tainty has  been  introduced.  That  tsj  does  tfie  samnlinq  start  before 
the  snectrum  reaches  the  highest  values  in  frequency?  If  the  siqnal 
had  been  the  phoneme  recorded  in  the  normal  voice,  this  would  not 
have  been  critical  since  the  snectrum  is  flat  for  a  reasonable  period 
of  time. 

Since  the  data  was  already  recorded,  it  was  not  possible  to 
record  on  an  adjacent  channel ,  any  form  of  siqnal  to  coincide  with 
the  start  of  the  sound  so  that  the  sample  process  could  have  been 
automatically  triqqered.  The  alternative  was  to  simultaneously 
monitor  the  signal  by  audio  and  visual  means  and  manually  triqqer  the 
sampling  process.  Obviously  this  technique  has  flaws  and  the  exact 
point  of  entry  into  the  data  stream  is  uncertain. 

Even  with  the  ambiquities,  the  fact  that  the  formants  are  shifted 
un  in  frequency  due  to  the  forced  sound  level  is  apparent.  The  first 
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FIG,  24 
Transient  Spectrum,  With  Muzzle,  0-3  KHz. 
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FIG.  25 

Transient  Spectrum,  80  feet,  0-3  KHz. 
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formant  can  be  estimated  at  840  hertz,  in  aqreement  with  Finure  15. 
The  second  formant  is  centered  at  1720  hertz  and  the  third  at  2550 
hertz.  Bandwidth  is  extremely  hard  to  estimate. 

2.  With  Oral  Cavity 

A  sample  spectrum  for  the  recordings  made  usinq  the  commun- 
icators with  the  oral  cavity  is  shown  in  Fiaure  20  (a-c).  The  litera- 
ture concerned  with  small,  non-radiatinq  enclosures,  such  as  the  oral 
cavity  used,  states  that  in  theory  the  effect  of  the  muzzle  is  to 
raise  the  formants  in  frequency;  however  in  reality,  the  distortion 
is  much  more  complicated.  "...  The  formants  may  become  indistinct 
rather  than  shift  in  any  direction..."  [21].  These  observations  are 
partially  borne  out  in  the  spectra  referenced  above.  F?  does  appear 
to  be  indistinct  as  does  F3.  F-,  contains  the  major  concentration  of 
power,  the  center  being  shifted  down  instead  of  un  in  frequency.  The 
inability  to  massage  the  data  to  make  concrete  observations  is 
reqretable. 

3.  Increased  Ambient  Pressure  With  Oral  Cavity 

A  sample  spectrum  for  the  recordings  made  at  80  feet  is  shown 
in  Figure  26  (a-c).  The  effect  of  increased  pressure  is  well  documented 
[1,5,6].  The  formants  should  shift  upward  in  freouency  in  a  non-linear 
fashion,  "...the  first  formant  is  affected  more  than  the  second 
formant.  It  is  also  observed  that  the  intensity  level  of  voiced 
sounds  increased  with  increasing  ambient  pressure  but  that  typical 
noise  sounds,  such  as  fricatives,  and  the  burst  part  of  unvoiced  stons 
display  a  drop  in  intensity  relative  to  voiced  sounds  ..."  [6]. 
This  referenced  source  cites  a  theoretical  expression  for  computinq 
the  upward  shift  in  frequency  of  the  first  formant. 
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FIG.  26(a) 
Phoneme  /AE/,  80  feet,  0-1  KHz 
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FIG.  26(b) 
Phoneme  /AE/,  80  feet,  1-2  KHz 
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FIG.  26(c) 
Phoneme  /AE/,  80  feet,  2-3  KHz. 
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F-  2  =  F  2  +  F  2  (p-1) 
In      11     wl   v   ' 

where  F,   is  the  frequency  of  the  first  formant  at  n  atmosnheres 
pressure,  F, ,   is  the  first  formant  at  one  atmosDhere,  and  F  ,  is  the 
resonance  frequency  of  the  closed  vocal  tract  at  one  atmosphere. 
AoDlvinn  this  expression  to  the  situation  under  discussion,  namely  a 
depth  of  80  feet,  the  theoretical  first  formant  shfft  can  be  calculated. 
Assuminq  a  resonant  frequency  of  170  hertz,  and  usinq  the  estimate  of 
840  hertz  for  F,,,  with  p  =  3.42  atmosnheres,  F,  is  comnuted  to  be 
378  hertz,  a  shift  of  38  hertz.  Examininq  the  spectrum  for  80  feet 
in  Fiqure  26,  the  first  formant  can  be  rouqhly  estimated  to  be  910 
hertz,  siqnificantly  larqer  than  expected.  The  source  of  this  shift 
is  uncertain;  it  could  be  related  to  the  increased  effort  of  the  diver, 
as  well  as  to  the  increase  in  pressure.  F~  and  F,  are  indistinct, 
the  power  distributed  in  many  frequencies  throuqhout  the  band. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  original  objective  of  isolating  the  distortion  of  the  under- 
water communicators  was  not  accomol ished.  Ths  problem  narrowed  to 
develonina  a  tool  to  examine  the  sneech  power  snectrum  in  detail  with 
high  resolution.  It  is  felt  for  non-noise  like  Dhonemes  that  the 
FFT  is  the  tool  available  to  the  analyst  to  perform  this  high  resolu- 
tion analysis.  Some  of  the  final  massaged  output  reinforces  this 
belief.  Line  snectra  can  be  the  final  result  of  speech  wave  forms  if 
a  number  of  the  troubles  encountered  are  corrected  in  subsequent 
research.  To  assist  in  future  work  a  summary  of  the  troubles  encoun- 
tered and  recommendations  for  correction  are  submitted. 

A.  DATA  COLLECTION 

1 .   Equipment 

The  proper  equipment  must  be  available  to  sunport  the  research 
in  all  phases  of  data  collection  and  remain  available  for  recall  during 
the  orocessinq  phase.  This  affects  a  number  of  topics,  ranqe  location, 
etc.,  but  is  best  illustrated  by  the  necessity  for  long  term  availabil- 
ity of  the  underwater  communicators.   In  this  experiment,  the  communi- 
cators were  on  short  term  loan  and  a  number  of  experimental  technigues 
were  hastily  contrived  to  support  the  data  collection. 

The  proner  divina  eguipment  must  be  made  available  to  suooort  the 
collection  effort.  Diver  comfort  in  cold  water  requires  the  addition 
of  heated  suits  to  prevent  changes  in  speech  due  to  loss  of  heat.  The 
orientation  problem,  discussed  later,  calls  for  a  helmet  or  hood  tvpe 
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lantern.  The  type  developed  for  Sea  Lab  III  by  Ben  Saltzer  of  Naval 
Undersea  Warfare  Center,  Pasadena,  California  is  recommended.  Short 
term  loan  of  nersonal  equipment  and  use  of  outdated  equiDment  compli- 
cates the  overall  effort  and  reduces  efficiency. 

2.  Tape  Recorder 

The  tane  recorder  deserves  special  mention;  the  choice  is 
critical.  The  PI  6200  was  chosen  because  it  could  be  adonted  for 
dc  power  and  nrovided  the  necessary  linear  frequency  translation  due 
to  the  slow  sample  rate.  The  PI  6200  tape  recorder  is  delicate  and 
required  special  handlinq,  but  is  suitable  for  use  at  a  fixed  recordinq 
site. 

3.  Range  Location 

The  ranqe  location  is  critical.  The  attempt  to  set  un  a  ranqe 
in  Monterey  Bay  durinq  the  worst  weather  in  years  cost  many  man  hours. 
The  ranqe  was  completed  and  used;  however  theft  by  local  fishermen 
negated  the  effort.   If  Monterey  Bav  is  to  be  used  in  the  future  (it 
should  not!)  the  summer  months  provide  the  best  ocean  conditions,  and 
notices  should  be  out  out  to  all  fishermen  askinq  them  to  cooperate  in 
avoidinq  destruction  of  the  ranqe.  The  Monterey  Harbormasters  Office 
will  assist  in  this  resoect.  Time  and  effort  can  be  saved  if  a  better 
location  is  used,  specifically,  Staqe  I  at  NSDRL,  Panama  City,  Florida. 

4.  Diver  and  Receiving  Transducer  Orientation 

Diver  orientation  must  be  closely  controlled  as  must  the 
orientation  of  the  receiving  transducer.  The  luna  cavities  and  wet 
suit  cause  the  transmission  of  the  communicators  to  be  hiqhlv  direc- 
tional.  If  the  diver  is  not  oriented  properly,  the  sianal  strength  will 
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be  reduced;  similiarly  the  WQC-1A  transducer  is  directional,  so  must 
be  oriented  to  nrovide  maximum  siqnal  reception.   In  Monterey  Bav  the 
divers  orientation  was  fixed  by  a  nlatform  allowinq  him  to  sit  down. 
No  attemnt  was  made  to  orient  the  transducer  in  the  ocean  ranqe.  At 
San  Vicente  Reservoir  the  qranite  walls  and  concrete  dam  face  required 
that  the  receiving  transducer  be  oriented  to  oreclude  reflections.  The 
diver  was  oriented  using  a  wrist  comnass  to  maintain  the  nroner  line  to 
the  receivinq  transducer.  The  diver  held  a  flashlight  in  one  hand  to 
assist  in  readinq  a  word  list  in  total  blackness,  while  he  clunq  to  a 
safety  line  to  control  his  depth  and  location.   It  is  no  wonder  the 
result  was  fadinq  due  to  diver  drift  from  the  nroner  diver/transducer 
line.  A  setun  of  the  DICORS  type  would  materially  assist  in  solvinq 
the  orientation  problem  [5]. 

5.  Safety 

Additional  divers,  preferably  with  a  technical  backciround,  are 
needed  for  safety  and  convenience.  Throuqhout  this  research,  the  divers 
violated  normal  safety  standards  by  diving  alone.  This  left  one  man 
at  the  surface  to  control  and  monitor  the  myriad  of  controls  for  re- 
cordinq.  The  ability  of  the  topside  man  to  monitor  the  data  collection 
with  confidence  was  hoDeless. 

6.  Spurious  Noises 

The  attenuation  of  noises  associated  with  tane  recorder  heads, 
interference  of  60  cycle  nature,  surface  and  bottom  reflections  must 
be  constantly  monitored.  Data  should  be  recorded  only  if  optimum 
conditions  exist. 

7.  Forced  Speech 

The  data  collection  must  suooort  the  data  orocessinq.  The  verv 
fact  that  two  indeoendent  projects  were  undertaken  reduced  the 
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effectiveness  of  either  project.  The  method  of  collectlnq  data  for 
intelligibility  tests  is  vastly  different  than  that  to  support  the 
speech  processing  -  the  subject  of  this  paper.  Intelligibility  tests 
are  optimized  when  the  diver  raises  the  level  of  his  voice  and  makes  a 
determined  effort  to  speak  clearly.  However,  as  mentioned  repeatably, 
the  task  of  isolating  the  problems  of  distortion  by  use  of  the  power 
spectrum  requires  the  diver  to  speak  in  the  same  manner  throughout. 
The  line  spectrum  desired  is  most  nearly  approached  when  the  voice 
is  articulated  in  a  normal  manner.  This  implies  that  instead  of 
forcina  the  voice  on  the  surface  to  agree  with  the  technique  used  for 
clarity  at  depth,  the  voice  must  be  used  in  a  normal  manner  at  each 
stage  of  the  experiment,  a  reverse  of  what  was  done. 

B.  DATA  PROCESSING 

The  processing  technique  must  be  thorouqhly  understood  and  be  used 
successfully  prior  to  the  collection  of  the  data.  Lack  of  a  hybrid 
interface  and  limited  availability  of  the  underwater  communicators 
precluded  these  conditions  being  met. 

1 .  Sample  Rate 

The  hybrid  computer  facility  should  modify  the  write  scheme 
to  ensure  that  the  A-D  conversion  rate  is  high  enough  to  accommodate 
speech  signals.  The  desired  rate  is  at  least  ten  kilohertz.  The 
analog  computer,  the  COMCOR  5000,  is  presently  rated  at  25  kilohertz; 
the  limiting  item  is  the  dialtal  computer,  the  SDS  9300.  The  limita- 
tion of  2.5  kilohertz  forced  the  conversion  of  the  data  at  a  playback 
speed  of  3.75  ips,  one-tenth  the  record  speed.  In  theory,  no  error 
is  introduced  but  it  is  desirable  to  process  at  the  record  speed. 
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2.  Filterinq 

Filterinq  of  the  siqnal  prior  to  A-D  conversion  should  be 
variable  and  of  hi qh  quality.  The  Krohn-hite  filter  is  suitable  for 
such  work  and  in  many  ways  is  superior  to  a  fixed  BandDass  filter 
usinq  active  devices.  A  voltaqe  controlled  bandpass  active  filter 
would  also  be  suitable. 

3.  Interfacing 

Converting  the  seven  track,  24  bit  tane  to  a  nine  track,  32 
bit  tane  was  a  particular  stumblinq  block,  requirinq  lonq  hours  to 
overcome.  Anpendix  B  is  included  to  assist  in  overcominq  this 
problem  in  the  future. 

4.  Output  Format 

The  output  from  the  IBM  360  can  be  displayed  in  a  variety  of 
ways;  however,  the  CalComp  Offline  plotter  is  most  suitable  for 
analysis.  A  more  efficient  draw  routine  would  make  processinq  easier 
and  save  compute  time  and  core  memory. 
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APPENDIX  A  GLOSSARY 

PHONEME—  one  of  a  finite  (40)  number  of  sounds  that  characterize 
speech.  Six  major  subqrouns  (classes)  are  shown  in  Table  VII. 

FRICATIVE—  a  noise-like  sound,  primarily  produced  by  blowinq  air  over 
the  teeth,  as  in  "s",  "sh\  and  "f.M 

STOP—  a  sound  produced  by  holdinq  the  breath  and  exolodinq  it  in  a 
burst  of  noise,  as  in  "p",  "t",  and  "k." 

VOWEL—  a  steady  state  sound  as  in  all  common  vowels. 

GLIDE—  a  sound  formed  by  continually  movinq  the  articulators 
(mouth,  tongue,  etc.),  as  in  "1",  "r",  and  "w." 

NASAL-  "m"  and  "n"-like  sounds. 

DIPTHONG—  two  or  more  vowels  sounded  toqether 

PLOSIVE-  a  stop 

FORMANT—  the  resonant  frequencies  of  the  vocal  tract. 

AFFRICATIVE—  a  nhoneme  beginning  with  a  stop  and  ending  with  a 
fricative,  as  "ch"  in  the  word  "church." 

PITCH—  a  qualitative  dimension  of  hearinq  related  to  the  hiqhness 
or  lowness  of  a  tone  and  the  sound  intensity  (with  increased 
intensity,  a  relative  chanqe  in  pitch  can  occur  with  a  chanqe 
in  frequency. ) 
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STOPS 

/p/  nat 
/t/  to 
A/  key 


/b/  be 
/d/  day 

/q/  qo 


TABLE  I  [21] 

FRICATIVES 

/h/ 

house 

hi 

vest 

/f/ 

fee 

1*1 

then 

/£/ 

thin 

/z/ 

zoo 

/s/ 

see 

1)1 

azure 

II 

she 

NASALS 

/m/  me 
/n/  need 
I'jl     song 


GLIDES 

/w/  we 

151  you 

/r/  read 

/V  "let 


VOUELS  (includina  dipthonqs) 


IM 

eve 

/AE/ 

at 

M 

boy 

/u/  foot 

HI 

it 

IV 

ask 

In/ 

not 

(Ml     boot 

le/ 

hate 

/Al 

up 

lol 

all 

/ai/  I 

Icl 

met 

leV 
10/ 

say 
father 

IV 

obey 

/au/  out 
/ou/  go 
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APPENDIX  B 

The  nurnose  of  the  anpendix  is  to  qive  a  brief  but  complete 
description  of  the  mechanics  involved  in  diqitizina  data  and  inter- 
facino  from  the  SDS9300  to  the  IBM  360/67. 

A.  ANALOG  TO  DIGITAL  CONVERSION 

The  analoa  sianal  is  the  input  to  the  Comcor  CI  5000.  For  all 
siqnals,  a  bandnass  or  lownass  filter  network  is  necessary  to  prevent 
aliasina,  a  major  problem.  The  upper  cutoff  frequency  should  be  the 
hiqhest  frequency  of  interest.  The  filtered  siqnal  is  pluqqed  into  an 
amplifier  input  on  the  CI  5000  and  amplified  by  the  hiqhest  possible 
factor  within  the  +100  volt  limits  of  the  operational  amplifiers. 
This  is  required  to  yield  the  hiqhest  possible  sianal  to  noise  ratio. 
The  output  of  the  amplifier  is  the  input  to  Trunkline  500(T500)  on  the 
analoa  patchboard.  This  links  the  analoq  to  the  diqital  computer. 

The  samplina  frequency  can  be  selected  up  to  and  includinq  2500 
hertz  for  the  present  diqital  maqnetic  tape  write  scheme.  If  the  analoq 
clock  pulse  is  used,  a  selection  of  10,  100,  or  1000  hertz  is  available 
on  the  logic  patchboard  and  should  be  connected  into  the  Master  Clock 
Input  (MC  IN).  The  output  of  the  Master  Clock  (MC)  is  pluqqed  into 
Trunkline  210  (T210)  on  the  logic  patchboard. 

If  a  variable  samolinn  frequency  is  desirable  or  if  a  rate  differ- 
ent than  what  is  available  on  the  loqic  board  is  required,  an  oscillator 
can  be  pluqqed  into  a  comoaritor  input,  (say  C03)  on  the  analog  patch- 
board. The  outnut  of  the  comoaritor  on  the  logic  patchboard  is  pluqqed 
into  Master  Clock  In  (MC  IN).  The  output  of  the  Master  Clock  (MC)  is 
connected  to  Trunkline  210,  (T210). 
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The  samplinq  proqram  for  the  SDS  9300  is  titled  SAMPL,  and  was 
developed  by  Lt.  Lynn  Dorrian  for  the  SDS  930  [14].   It  was  modified 
for  the  SDS  9300  by  Lt.  Jerry  Post  and  Mr.  R.  Limes  [15].  The  pronram 
will  sample  one  analoq  sianal.  The  data  is  stored  in  buffers  and 
written  on  maanetic  tape  from  these  buffers.  It  is  desinned  to  fill 
one  buffer  while  the  other  buffer  is  dumninn  onto  the  mannetic  tane. 
Since  the  write  scheme  for  the  magnetic  tane  is  inherently  slow,  the 
samnlinq  rate  is  presently  limited  to  a  maximum  of  2500  Hz.  If  only 
two  buffers  of  data  is  desirable,  this  rate  can  be  slightly  exceeded. 

The  procedure  is  to  call  the  Subroutine  SAMPLE (NREC,NSAMP, IBUFO, 

IBUF1,NTAPE,IT0G)  where: 

NREC  =  number  of  records/run. 

NSAMP  =  number  of  samples/record. 

IBUFO  =  first  data  buffer. 

IBUF1  =  second  data  buffer 

NTAPE  =  maqnetic  tape  unit  number 

ITOG  =  0  if  IBUFO  is  currently  beinq  filled. 
=  1  if  IBUF1  is  currently  beinq  filled. 

After  the  proqram  has  been  accepted  by  the  SDS  9300  and  the  parameters 
initialized,  the  diqital  computer  waits  for  the  analoq  computer  to  ao 
to  "compute."  The  compute  switch  can  be  manually  trioqered  or  trionered 
by  a  control  siqnal . 

Once  the  A/D  routine  commences,  it  continues  until  either  the  analoq 
qoes  out  of  "compute"  or  the  number  of  records  of  data  has  been  taken. 
Insure  after  the  proper  number  of  records  have  been  taken  that  the  key- 
board is  out  back  into  reset,  if  triqqerina  manually.  Control  returns 
to  the  dioital  proqram.  The  diqital  proaram  waits  for  instructions. 
If  Sense  Switch  6  is  liqhted,  the  proqram  returns  to  initialization  and 
waits  for  the  analoq  to  aoain  qo  to  compute  to  diqitize  a  new  run  of 
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data.  If  Sense  Switch  5  is  li anted,  the  maqnetic  tane  is  rewound  to 
load  noint.  If  Sense  Switch  4  is  liahted,  the  End  of  File  condition  is 
entered  on  the  tane.  Sense  Switch  3  will  cause  the  line  nrinter  to 
nrint  the  last  2  records  of  data  in  a  nrint  plot  routine. 

Enclosure  (1)  is  a  diaqram  of  the  CI  5000  analoq  and  Ionic  natch- 
board  configuration  for  the  A/D  conversion.  Enclosure  (2)  is  the 
nronram  listinp  for  the  Fortran  IV  routine  utilized. 

Below  are  listed  a  number  of  hints  to  save  effort  and  time. 

1)  Use  a  maanetic  tane  provided  by  the  IBM  360/67  comnuter  facility 
as  it  has  a  nermanently  marked  load  point  with  the  correct  snacinq. 
The  tape  unit  number  is  1  and  should  be  set  on  "automatic." 

2)  Ensure  that  dipital  proqram  uses  fresh  cards,  narticularlv  the 
RTM  BOOT  card;  a  feed  check  on  the  card  reader  is  often  sensed  as  a 
Droqram  error. 

3)  Follow  basic  9300  instructions.  If  the  control  console  fails  to 
start  oroaram  feed,  turn  off  the  system  and  start  over.  This  "immediate 
action  battle  drill"  could  save  manv  hours. 

4)  Use  multichannel  oscilloscope  on  the  analoq  console  to  monitor  the 
ore-filtered  sianal,  the  post-filtered  sianal  and  the  amnlified  sianal. 
This  will  allow  a  visual  check  to  ensure  proner  analon  action. 

5)  Monitor  hybrid  interface  console  durinq  diqitizinn  to  ensure 
that  the  system  is  workinq. 

6)  Make  at  least  the  first  two  and  last  two  runs  from  a  known  analoq 
sinnal  (known  frequency  and  amplitude).  This  is  necessarv  to  ensure  the 
conversion  to  the  IBM  360  is  complete  and  correct.  An  additional  suq- 
nestion  is  to  insert  a  run  of  known  data  periodically  throuqhout  the 
data  of  interest  to  also  allow  a  check  of  the  conversion. 

7)  Read  NPGS  Comnuter  Facility  Technical  Mote  0211-64  prior  to 
attemntinn  A/D  conversion. 
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B.   CONVERSION  TO  THE  IBM  360/67 

The  proqram  used  to  convert  from  seven  track,  24  bits  to  nine  track, 
32  bit  data  is  listed  in  Fnclosure  (3).  The  proaram  is  both  Fortran  IV 
and  assembly  lanquaae.  The  control  cards  are  listed  below. 
Card  No. 

1  //Name   (standard  job  card)  ... 

2  //FOURIER  EXEC  FORTCALG, TIME. G0=05, REGION  .GO=160K 

3  //FORT.SYSIN  DD  * 

(main  FORTRAN  IV  PROGRAM) 
/* 

4  //ASM.SYSIM  DD  * 

(ASSEMBLY  LANGUAGE  PROGRAM) 

5  //G0.FT02F001  DD  UNIT=0C0,VOL=SER=NPSYYY,LABEL=( ,BLP) ,      * 

6  //  DISP=(OLD,KEEP),DCB=(DEN=0,RECFM=F,BLKSIZE= 
4096) 

7  //G0.FT04F001  DD  IINIT=2400,V0L=SER=NPSXXX,LABEL=(  ,SL)  ,DISP= 
(NEW, KEEP),      * 

8  //  DCB=(DEN=2,RECFM=F,BLKSIZE=4096),DSNAME= 


Card  #1,2,3,  are  the  standard  control  cards  for  the  NPGS  Computer 
facility.  Card  4  indicates  the  followinq  proqram  is  in  assembly  lan- 
quaqe.  Card  5  and  6  tell  the  computer  where  the  seven  track  data  is 
located  and  its  format.  For  example,  the  card  above  says  the  7-track 
data  is  located  on  tape  NPSYYY,  to  be  placed  on  tape  drive  02,  the 
label  process inq  is  to  be  bypassed,  with  fixed  format,  with  blocksize 
equal  to  4096  bits,  (1024  words).  Cards  7  and  8  tell  the  computer 
where  the  data  is  to  be  transferred.  The  data  is  to  be  written  on 
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tape  NPSXXX,  tape  drive  04,  with  standard  label,  fixed  format  in  blocks 
of  4096  bits  (1024  words).  The  1024  words  corresponds  to  the  number  of 
samples/record  in  the  Subroutine  SAMPL  on  the  SDS  9300  when  NSAMP=1024. 

Anticipate  difficulty  in  the  conversion  of  the  data.  The  tane 
drives  are  the  weakest  link  and  will  require  persistence  to  be  success- 
ful. The  first  record  of  the  data  will  probably  be  missed.  The  known 
data  will  allow  a  quick  check  of  the  system. 
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Enclosure(2)  to  Appendix  B 

c 

C  THIS  PROGRAM  SAMPLES 

C  RHNO  AND  HYATT    DIVERS  SPEECH 

C 

C 

r  *  TN  COLUMN  6  INDICATES  CONTINUATION  FROM  PREVIOUS  CA&D 

C  FOR  THIS  PROGRAM  LIST  ONLY 

c 

Apirtpan   GO.LS 

DTMFNSIPN    IBUFK1024) ,IBUF2(1024),BUF(2048,2),JXY<2), Y 
*M 513) 

NAMELTST    NCBLKS, LENGTH, ITAPE,ITOG 

ITA£>E  =  1 

NOPLK*=2 

LFNGTH=1024 

ITOG=0 

jxy(  n=i 

JXY(2)=2 

1  OUTPUT(lOl)    NOBLKS, LENGTH, ITAPE,ITOG 
WRITE( 102,104) 

104       FHRMATft    TYPF    IN    CHANGES. %/%    FOLLOW    WITH    AN    *    AND    RETU 
*RN. %) 
INPUT! 1C1) 
C 

101  FQRMAT(2I2) 

2  IF( ITAPE)3.4,3 

3  WRTTE( 1C2.102)ITAPE 

102  c°PMAT(f   CHECK    ON    TAPEM1S    ANO    CLEAR    PAUSE.*) 
PAUSE     10 

4  CALL    SAMPLf NCBLKS, LENGTH, IBUF1, IBUF2,ITAPE, ITOG) 

5  !F( ITOG.EO.O)GO    TO    5 
*            CONTINUE 

NT=0 

IF(  ITOG.EO.-DGO   TO    10 
7  TM  ITOG.EO.DGO    TO   7 

B  CONTINUE 

NT=1 

IF(  ITOG.EO.-DGO    TO    10 

Gn    TO    5 

10  WRITF(102.103) 

103  FOPMATf //5X4END    RUN...$MSET    DESIRED    SENSE    SWITCH. i) 
20         PAUSF    20 

IF( SENSE    SWITCH    3)11,12 

12  IF(  SENSE    SWITCH    4)13,14 
1*         I»=<  SENSE    SWITCH    5)    15.16 
16         IFf  SENSE    SWITCH    6)1,20 

11  ML=( 1-NT)*LENGTH 
ML=NT«LENGTH 

^0    IP    1=1, LENGTH 
BUFU  .2)  =  f 

BUF( ( f+LENGTH) ,2)=I+LENGTH 
BUF( I*ML,1)=IBUF1(I> 
IP  RtJFM  I+NL)  ,1)  =    IBUF2U) 

L4=LENGTH/2 
L=LENGTH*2 

CALL    L0NGPL0T(BUF,JXY,L,2048,1,0,0.0,0.0,0.0,0.0,L4,YS 

GO    TO    2C 

13  ENOFILE    ITAPE 
GO    TO    20 

1<5         REWIND    ITAPF 

GO    T0    20 

EMO 
Amcta<5300       SI,LO,GC 
t^AMPL    07.F 
* 

*  THIS    PROGRAM    SAMPLES    ONE    ANALOG    SIGNAL    ON    TRUNKLINE    0500 

*  A""    A    RATE    DETERMINED    BY    A    CLOCK    PULSE    APPEARING    ON   TRUNKLI 

*  0210.       ANALOG    TO    DIGITAL    CONVERSION    IS    PERFORMED    AND    THE 

*  0AT*     IS    STOPED    IN    BUFFERS    SET    UP    IN    THE    CALLING    PROGRAM. 

*  CALLING    SEQUENCE    IS... 
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*  CALL     <AMPL(NREC, NSAMP, IBURO, IBUF1 ,NTAPE, ITOG) 

*  WHF  P  p  . 

*  NRFC=NUMPER    CF    RECORDS/RUN 

*  N^MP  =  NU*BER    OF     SAMPLES  /RECORD 

*  IBMpO=PIRST    PATA    BUFFER 

*  iBUFi  =  <;FrnND  data   buffcr 

*  IBiJFO  =  FIRST    DATA    BUFFER 

*  BUFFERS    APE    FILLED    WITH    NSAMP    A/O    SAMPLES    ALTEPNA 
♦  TftY 

NTAPE  =  M^G    TAPE    UNIT    NUMBER 

IF=C,    NC    TAPE    IS    WRITTEN 
ITOG=0    IF     IBUFO    CURPFNTLY    BEING    FILLEO 
=1     IF    IBUF1    CURRFNTLY    BEING    ^ILLED 


* 
* 
* 
* 

* 

* 


* 
TTG 


AFTER  INITIALLING,  THE  SUBROUTINE  WAITS  FOR  THF  ANALO 
*G  TO 

GO  TO  COMPUTE.  END  OF  RUN  TS  CAUSED  PY  ANALOG  GOIN 
*G  OUT 

OF  COMPUTE  OR  WHEN  NPEC  RECOROS  OF  DATA  HAVE  B«=rN  TAKE 
*N. 

IN  EITHFR  CASE,  ITOG  IS  SET  TO  -1  TO  NOTIFY  THF  CALL  IN 
*G 

ROUTINE  OF  THE  ENDRUN  CONDITION. 

NOTF  THAT  AFTER  INITIALIZING  THE  A/D  SAMPLING,  CONTROL 
*IS 

RFTURNED  TO  THE  CALLING  ROUTINE  AND  INTERRUPTS  CALL  TH 
*E 

A/D  INPUT  ROUTINE. 


NPEC 
NSAMP 

BUP<=0 
BUFFI 

TAOC 

TOGGL 

* 

a 

P 

XI 
X2 
X3 


OPO 
BRM 
R7E 
ITG 
!TG 

TTG 

ITG 
TTG 
ITG 
*G  FILLED 

EOU 

EOU 

EOU 

EOU 

EOU 

EOU 

EOM 

POT 

LDA 

STA 

LDA 

STA 

LDA 

ETR 

SKU 

STA 

COPY 

MPG 

STA 

COPY 

MRG 

*TA 

COPY 

MPG 

STA 

COPY 

MPG 

STA 

LDA 

COPY 

SU« 

MPG 


0100000 
9SETUPN 


5 

-5 

4 

1 

2 

3 

033001 

=0 

♦TOGGL 

TIME 

=-1 

NOTAPE 

♦  TAPE 

=  0T 

=  0 

NCTA PE 

<A,B) 

WTBX 

WTBN 

IB, A) 

EFTX 

EFTN 

(B,A) 

TRTX 

TRTK 

(B.A) 

BTTX 

BTtn 

♦NSAMP 
(A,B) 
=  1 
PLACE 


NUMBER  OF  RECOROS/RUN 

NUMBER  OF  SAMPLES/RECORD 

FIRST  DATA  BUFFER 

SECOND  DATA  BUFFFR 

MAG  TAPE  UNIT  NUMBER 

TOGGLE. ..INDICATES  BUFFER  BEIN 

OR  ENDRUN  CONDITION. 


SETLINES  1  AND  2  FALSE 
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STA 

sna  tat 

c;ta 

C  UR  X2 

ropy 

(0VA) 

LLT 

14 

str 

OAT 

vpg 

*LC 

CT» 

ALCC 

STA 

A  1  C  1 

CPPY 

(0,M 

LOR 

PUcPO 

Lt  SP 

Q 

LLrP 

1 

\  R  SR 

10 

LI  SA 

«; 

"RG 

«LC0 

CJA 

ALCG 

rnoY 

(B.M 

PAT 

mot, 

STA 

CWO 

f  no  y 

(0,*) 
BUFFI 

LPB 

LLSP 

Q 

LL^P 

1 

LP  SB 

10 

Lt  rA 

5 

MPT, 

A  LCI 

STA 

A  LCI 

r*OPY 

(Bf  A) 

wr, 

P£T 

*TA 

CW1 

LP* 

RUFFO 

«=tr 

=0T777T 

*pp 

=  -1 

STA 

ORfr,o 

STA 

CO*1* 

lpa 

BUFFI 

«=TR 

=077777 

Ann 

=-1 

sta 

0»IG1 

LPA 

OCW 

mog 

=CP*M 

CTA 

svcv 

STA 

cw 

LP& 

BRMflO 

XMA 

040 

STA 

SV040 

LPA 

rbmsam 

XMft 

051 

STA 

SV051 

SKN 

NQTfipc 

BPU 

AG*  IN 

LPA 

ROMCL 

XMA 

010 

STA 

<V010 

LPA 

BRMFL 

XMA 

Oil 

STA 

SV011 

CXU 

TR^N 

CAT 

0 

RRU 

f-2 

FXU 

BTTN 

RRU 

t*2 

BPU 

AGS.  IN 

COM 

♦014000 

DOT 

^DACF 

rxu 

BFTN 

aga  im 

ST7 

♦TOGGL 

ST7 

TOGGT 

LPA 

=-1 

STA 

APFL 

LPA 

♦  NOEC 

APP 

=-1 
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*  FND    0 

*  SlIRPP 

* 

cT/ypT 


LOAPl 


RTCCT1 


PT^ST2 


STA 

PPM 
PRT 
SK«: 
RPU 
POM 
DOT 

PTR 

EOM 
BPR 

*=     TNI 

IJTTWP 

PZE 
STA 
STX 
c*N 
BPIJ 

OIR 
lda 

S*G 
BPU 
MPQ 
S*S 
RPM 

ST7 

PPM 

PPT 

FIP 

LDX 

RRX 

S*N 

RPU 

LDA 

$ta 

LHA 

STA 

ST7 

LDX 

SKN 

BPU 

EXU 

CAT 

RPU 

LDA 

ADM 

ADD 

SKL 

COPY 

STA 

EXU 

EXU 

POT 

BPU 


CCUNT 

C33001 

=0^0000000 

030010      ANALOG 

*-l 

C33C01 

=  0 


SETLINF    1    TPUE 
IN    COMPUTE    TEST 


033004 
SAMPL 
TIAL    SET    UP 
START    ENTFPED 


SVA 

SVX2,X2 
ADFL 
«-l 


FALSF 


ENABLE 
RETURN 


PATCHBOARD    INTERRUPTS 
TO    MAIN    PROGRAM 


ON    INTERRUPT    051 


♦TOGGL 

TOGGT 

LOADO 

COMM 

030010 
END  PUN 
AOFL 
C34C00 
CW 


THIS    PROCEDURE    TESTS 
TO    DETEPMINE    WHICH 
BUFFER    TO    LOAD 

ANALOG    IN   COMPUTF 

NO 


,X2 
»X2 


CURX2 

PTDLE 

4DFL 

t-1 

ORIGO 

COMM 

svcw 

c  w 

♦TOGGL 

SDATAT,X2 

NOTAPE 

INCR 

TRTN 

0 

t-2 

TIME 

*BUFF1 

=  1 

=  1000 

CO, 5) 

TIME 

WTBN 

A  LCI 

CW1 

I  NCR 


*       THIS    DIVIDES    THE    SUBROUTINE    INTO    BUFFERS 

LOADO 


RTEST3 


MPO 
<:ks 
BRM 
ST7 
EOM 
POT 
FIP 
LDX 
BRX 
SKN 
BPU 
LDA 
*TA 
LDA 


COMM 

030010 

ENDPUN 

ADPL 

C34000 

CW 

CURX2tX2 
RIDLE.X2 
ADFL 
$-1 

0RIG1 

COMM 

SVCW 


ANALOG 
NO 


IN    COMPUTE 
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CT/i 

r  v^ 

mtp 

•TOGGL 

L^X 

<:r»&TATtX2 

p  tc<;t4 

*KN 

NOT?PF 

ROM 

IMCR 

SXl) 

T«?TN 

r/»  t 

0 

ROIJ 

t-2 

L^A 

▼T  *F 

inn 

♦  RtJFPO 

Ann 

=  1 

?XL 

*1000 

rnpY 

(0,M 

CT* 

TJ  MC 

PXU 

WTRM 

*XU 

ALCO 

DOT 

CWO 

TMPC 

*KP 

CPJJNT         RL 

RP  11 

RTPLE          IP 

*r 

TP 

RDM 

PIN 

C    Tr>l   C 

ma 

CV* 

*TX 

flJPX2tX2 

LPX 

«VX2»X2 

ROC 

c  IM 

^  jo 

033000 

COM 

033001 

DOT 

^020000000 

<-*c 

C31001 

RPII 

♦  -1 

POM 

C33001 

nnr 

»o 

<:kp 

♦TOGGL 

"Pll 

t-1 

CAT 

c 

"PH 

*-l 

LP* 

rvo^o 

<T£ 

C40 

LHA 

CV0=51 

?T* 

051 

lpa 

CV010 

CTfi 

010 

LPA 

svon 

CT^ 

Oil 

PJC 

rpii 

pyoiF 

e"»n«»l|W 

HTC 

LP' 

r  pmh/ 

Crp 

i£77777 

CTi 

cc^m 

LP  A 

*-i 

I  nx 

riJRX2»X? 

rTj 

*C0MM 

Mnq 

CCM" 

RPX 

f-2,X2 

<"T7 

COUNT 

!  P.* 

rNoPijM 

fk'", 

■L0ADC 

RPH 

CT^T? 

£ 

RPI) 

RTP  cja 

*           r\i 

D    OF    « 

L'PCCUT!NF 

R3  '/I  PL 

RO  M 

PLUG 

•3  1   ijr. 

C/c 

ncr 

*»UH3 

svoic 

P7r 

«  V  *  1 1 

OJZ 

* 

r  ~^j 

PDPV 

10,14 

'">\|T 

POP  V 

)CK    CHIINT    I<     OEPHCEn. 
NPGATIVC,     ALL    PATA    HA?     BFF.M 
TPAN«FF.PRPp. 


CPTLINP    2    triic 

WAIT  POP  ANALOG 
^N  TEST  i. 
SPTLINP  2  FALSF 
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SPACE 

CONT 

150.0 
C100000 

new 

DAT* 

SVCW 

P?F 

cw 

P^E 

COMM 

PZC 
P7E 

SDATAT 

RES 

1 

OR  TOO 

RFS 

1 

OR  TGI 

RCS 

1 

DAT 

PES 

1 

ALC 

POM 

C14000 

ALCO 

P7E 

ALC1 

P7E 

rwo 

RFS 

1 

rwi 

RPS 

1 

TOGCT 

RES 

1 

COUNT 

RES 

1 

SVA 

P7F 

SVX2 

P7E 

CURX2 

P7E 

SVOAO 

RES 

1 

SV051 

RES 

1 

PLACE 

DATA 

077700000 

tjmc 

PZE 

NOT/>PE 

P7E 

BTTX 

RTT 

0,0 

B^N 

P7F 

WTBX 

WTR 

♦0,0,* 

WTRN 

°ZE 

EFTX 

EFT 

0,0,4 

EF^N 

D7F 

TRTX 

TRT 

0,0 

TR^N 

P7E 

BRM^AM 

BRM 

START 

BPMAD 

BR* 

A  OF  NO 

ADENO 

P7E 

SKP 

AOFL 

BPU 

t-l 

BPC 

*ADEND 

AD^L 

P7E 

CND 

4E0F 

ALO^D 

XM,MAP. 

6 

DATA 
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Enclosure(3)  To  Appendix  B 
THI*"    DpriGRAM    INTERFACES    THE    SDS    9300    TO    TH~     TBM360 


'/my^000c3    J^p    (0053.47FP) • «MYATT  • ,  MSGLFVaM ,  CL  &SS=G 

//-PUPTcr    FXT    ERTCAIGD,TI  *E. OP=10 , REGION. 00= 160 K 
//pppt.<:y«in  p^    * 

oimpn^icn  indata(102m  ,0ata(1024) 

r^rTPq=1Po#o/f 2**31-1 ) 
crwiMn    ? 

nn    j  =  lfK  WHERE    K    IS    THP    NUMBER    OF    BLOCKS    OF    102'+    EAC 

^°    31    J=l,133 

RcAn( 2,3,END=40,FPR=C0)     INDATA 
3    cnrMAT(  16(  64A4)  I 

CM  I     FP0M(  TNPATA) 

W°  T  TF  (  *  ,  7  o )     J 
70    copm'T    (  •  1 ',10Xf 'RECP9P    NO.=',I6) 

OP    I    1=1.102* 
I    ">A'A f I )=TNnflTA( [)*FACTOR 

WRITc(6,66)     (D«TA(I)  ,  1=1  ,1024) 
66    PnpMATf  1X,10E1P.2) 

WQTTr     fA?3)    pATA 

GP    to     31 
^C    WRT^F     (6.51)     J 

*1    fppmat    (  •0S5X, 'READ    EPPOR,     RECORD    N0.=«,I6) 
31    CPNTTNUE 

STOP 
^0    WRIT-     (ff4i)     J 
41    FPPm,»t    (i0S5XtfFNn    pc    TAPEt    RECORD    NO.  =  «,IM 

^TOP 

//-"^M.^Y^IN       DP        * 

=  PRM  fT|\pj     o 

*  SUBROUTINE    FORM(INDATA) 

St 

H  *    TMPATA    re    an    CRRAY    LENGTH    S°ECTFIED    BY    THE     IN^EX    V£LU 

*  THIS    SUBROUTINE    WILL    CONVERT    24    BIT    BINARY    WORDS    STORED    IN 

*  *«=    Tp     32    PIT    BINARY    WORDS    AND    PLACE    tHESE    SAM*    WnP0S     BA 

*  tC*     INTO    INDATA 
H 

fTM  1^,12,12(13)  THIr    SUBROUTINE    CONVERTS 

BALP  6,0  24    BIT    BINARY    WORDS    TO 

USING  *,6  32    BIT    WORDS 

UMNG  DATA, 7 

SP  7,7 

L  11,=F'1C24«  THIS    IS    THE    IND*X    VALUE 

L  12,0(1) 

O^P  L  2,NUM(12) 

1°  3,7 

S°DL  2,6 

coL  2,2 

rPDI  2,6 

^L  2,2 

^°DL  2,6 

"L  2,2 

:rdl     2,6 

CT  3,NUM(12) 

L*  12,M12) 

BCT  ll.LOPR 

LM  2,12,28(13) 


MVI  12(13), X'FF* 

RCR  15,14 

>AT/  PSCCT 

ENO 
^/gc.fto2pooi  pd  unit*oco,vol=ser*divers  ,L^3EL=(,BLP)  , 
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//  DISP»(OLD,KFEP)  ,DCB«(DEN«0,«ECFM=FfBLK<;  IZF=40 

//nn.FTQAFOOl    DO    UNIT=2400,VOL«SER=NP$234fl_ABFl  =  (,SL),  01$°=  ( 

iNFW,KEFP),    * 
//  0CB=(nEN=2,RECFM»FfflLKSIZF=40°6),0SNAME=HYA00 

$053 
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COMPUTER  PROGRAM  HARM 

PROGRAM  COMPUTES  PGWER  SPECTRUM  FOR  8192  DATA  POINTS 
C 

%     IN  COLUMN  6  INDICATES  CONTINUATION  FROM  PREVIOUS  CARO 
C     FCP  THIS  PROGRAM  LIST  CNLY 
C 

//MYA20C53  JOB  ( 0053 ,47FP  )  ,  •  MYATT      '.MSGLEVEL-l.CLASS-D 
//FOURR  EXEC  FORTCLGP,PARM.F0RT*»LINECNT*75,,TIME.G0*02,REGI 

fON.GO=300K 
//FORT.SYSIN  DO  * 

REAL*4LABEL( 10)/40H0NE  TWO  THREFOURFIVE  S  IXSE VNE IGTNI N 
*E  TEN       / 

RFAL*8ITITLE(12)/96H      POWER  SPECTRUM      J.M.  MYATT 
1    M**2   111-118  W    MUZZLE  SINGING 

DIMENSION  XG(42C) ,D< 1260 ).Y < 8200 ).B ( 1260 ) 

DIMENSION  S(8200) , I NV ( 8206) . M(3 ) ,DATA( 1024) 

COMPLEXES  A<8192,i,i),Z< 1260,1,1 ) 

CALL  CANCELI2) 
C 
C 

ITYPE=2 

KK  =  13 

TIME  =ITIME(0)*0.01 

N1=2**KK 

NPT=N1 

NY  =  1 

NN'=420 

NZT=420 
C  PROGRAM  CCMPUTES  8192  COEF. 

C       PLOT  GOES  TO  3000  HERTZ 

IMAX=1 

BMAX=0.0 

KA=2 

KC  =  2 

KC  =  2 
C       SAMPLING  PATE  IS  2048  CPS. 

XN  =  N1 

DT=1. 0/2048.0 

BW=  1.0/(2.0*DT) 

DELTF=  1.0/(XN*DT) 

N2=l 

N3=l 

M(1)=KK 

M<2)=0 

M(3)=0 

KB=0 
C  NN=  RUNNING  COUNT  OF  POINTS  BEING  PLOTTED 

C  NZT=  NUMBER  OF  POINTS  PER  GRAPH 

C  DELTF  =  FREQUENCY  RESOLUTION 

C 
C 
C        J  IS  f  OF  BLCCK  OF  DATA  ON  MAG  TAPE  BEING  READ 


C 


DC  700  J=1,L18 

I0=KB*1024 

READ  (4,50)  (DATAU  ),  1*1,  1024) 
50  FORMAT(  U(64A4)) 

IF(J.LT.111)G0  TO  700 

DO  777  1=1,1024 
777  Y( IO-*I)=DATA<I  ) 

KB=KB+1 
70C  CONTINUE 

DO  2  11=1, Nl 

DO  2  12=1, N2 

DO  2  I3=1.N3 
2  A<  U,I2,I3)=Y(I1) 

WPITE(6,40)DT, J.Nl.NZT 
40  F0RMAT(////3X, 'SAMPLING  INTERVAL*  •  ,F10.6  ,  10X, •  SOUND 
i/AE/       ■ 
*,10X,«  BLOCK  =  »,I6,5X, 'DATA  POINTS • ,217/  ) 

WRITE(6,66)BW, DELTF 
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66  F0RMATC//6X, •BANDWIDTH  =  • ,  F15.  6,  10X,  •  DELTA  F  »  »,F10. 
$6) 

WRITE(6,10)  (  ((  II,        A(Il,I2,I3),I3=l,N3),I2  =  l,N2),  I 
$1=1, Nl) 

10  FORMAT<//»  AMPLITUDE  CF  FOURIER  FUNCTION  TO  BE  TRANSFO 
$RMED»/(4(1X, 

♦  •AC,  15,      • )  =  «,1X,2F7#3,4X))  ) 

CLOCK=ITIME(0)*C.6l 

CALL  HARM(A,M,INV,S,-1,IFERR) 

DC  34  I1=1,NZT 
34  XG(I1)=I1-1 
C    KO  IS  TEST  FOR  SMOOTHING 

IF(KO.LT.l)GO  TO  800 

Z(1,1.1)=A<1,1.1) 

B<l)*CABS(Ztl.l.l)) 

DO  801  11=2,1260 

12*11-1 

13=11*1 

ZU1.1,  l)=-0.25*(A<I2,l,l)*AU3,i,l))*0.5*A(Il,l,l) 

B(Il)«CABS(Z(tl, 1,1)1 

B(I1)=B(I1)**2 
801  CONTINUE 

WRITE<6.965) 
965  FORMAT! 1H0,»   FIRST  SMOOTHING  FUNCTION  USEO  •) 

GO  TO  820 
800  DO  33  11=1,1260 

B(  I1)=CABSU(I  1,1,1)  ) 
33  B(I1)  =  BU1)**2 
820  DO  11  11=1.420 

11  XG(I1)=XG(I1)«DELTF 
Jl=  0 

C      KC  IS  SECOND  SMOOTHING  FUNCTION  CHECK 
C 

IFCKC.GT.UGO  TC  987 
990  DO  666  11=1,1260 

IF<B(I1).LT.BMAX)G0  TC  666 

BMAX=B(I1) 

IMAX=I1 

666  CONTINUE 
WRITE(6,667)BMAX,IMAX 

667  FORMAT(1HO,»  8MAX    =    •,F7«3,«       IMAX    *    SI6I 
DC    78    J3=l,3 

J2=    1 
C  KA    IS    CHECK    FOR    LOG10    PLOT 

IFCKA.LT.DGO    TO    840 
DO    850    I1=NY,NN 
B(I1)=B(I1I/BMAX 
IF(B(I1).LE«0.002)GO    TO    1 
D(J2)  =  10.0*AL0G10(BUi)) 
IF(D(J2UGE.-27.0)G0    TO    850 
1    D(J21=-27.0 

850  J2=J2*1 

851  CALL  DRAW(K'ZT,D,XG,0,0,LABEL(J3),ITITLE, 3.  ,7.0,0,0,0,0 
$,9,15,1.LAST) 

GO  TO  855 
840  DO  77  U=NY,NN 
D(J2)=B(I1) 
J2=J2+1 

77  CONTINUE 

CALL  DRAW< NZT,D,XG,0,0, LABEL! J3), IT ITLE,0 ,7. 0,0, 0,0,0, 
$9, 15,1, LAST) 
C      NOTE$$t$$$$  X  AND  Y  AXIS  ARE  REVERSED 
C 

855  WRITE<6,20)(<( II,       A ( II ,1 2, 13) , I3«l ,N3 ) , 12*1 ,N2 ) , I 
$1=NY,NN) 
2C  FORMAT(//fFCURIER  TRANSFORM  A ( 1 1, 12, 13 ) • / (4( IX, • A ( •  15, 
$      »)*»,1X 
*,   2F7.3,4X))) 
NY=NN+1 
NN=NN+420 

78  CONTINUE 
ITYPE=2 
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GO  TO  35 
987  D(1)=B(1» 

TC  988  11=2, 1259 
12=11-1 
13=11*1 

DCI1I"  0«5*B(UI*0«25*tBf  121+Bf  1311 
98P  CONTINUE 

DC  989  11=1,1259 
B(  I1)=D( III 
989  CONTINUE 

WFITE(6.966I 
966  FORMAT(lHO,»    SECOND  SMOOTHING  FUNCTION  USEDM 
GO  TO  990 
35  CONTINUE 
END 
//GC.FTCAF001  DD  DSNAME  =  ^YA0053,UNIT=*2400,  DCB*  (RECFH«F  ,  BLKS  I 

tZE=4096),   * 
//  V0LUME=(PRIVATE,SER-NPS234I 
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